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I.  INTRODUCTION 

The  purpose  of  this  report  is  to  update  previously  reported  re- 
search progress  [l]  on  the  identification  of  naval  vessels  using  sub- 
structure complex  natural  resonances.  A multitude  of  details  on  a harmonic 
frequency  reflectivity  range,  complex  natural  resonances  of  structures, 
valid  electromagnetic  models  of  naval  vessels,  valid  electromagnetic 
models  of  the  sea  surface  and  a "matched-filter"  type  interrogating  radar 
signal  have  been  given  in  previous  reports  [1,2]  and  will  not  be  repeated 
here.  One  should  note,  however,  that  the  use  of  a conducting  ground 
plane  to  model  the  sea  restricts  the  frequencies  that  can  be  simulated 
to  some  highest  full  scale  frequency  in  the  neighborhood  of  25.0  MHz. 

The  model  scale  factors  used  (500:1  or  700:1)  allow  good  geometrical 
and  electromagnetic  simulation  but  it  is  not  possible  with  the  present 
system  to  obtain  valid  scattering  data  at  frequencies  higher  than  25.0 
MHz.  In  fact,  if  the  sea  surface  must  be  considered  rough  it  is  doubt- 
ful that  any  physical  ship-sea  model  can  be  used.  This  problem  does 
not  limit  the  feasibility  studies  of  the  present  effort  but  does  dictate 
that  our  methods,  if  successful,  cannot  be  conveniently  tested  at  fre- 
quencies greater  than  25.0  MHz.  The  present  feeling  is  that  identification 
at  frequencies  much  greater  than  those  presently  being  tested  would  have 
to  be  based  on  certain  substructure  complex  natural  resonances  whose 
extraction  would  be  seriously  affected  by  the  sea  state  and  wind  state 
as  well  as  other  variable  factors.  The  resonance  frequencies  of  interest 
are  those  with  the  vessel  in  the  presence  of  the  sea.  In  this  case  the 
quasi- invariance  of  the  resonances  might  be  lost  and  identification 
seriously  distorted.  The  interested  reader  is  referred  to  [1]  for 
pertinent  details. 

It  has  been  demonstrated  [1]  that  complex  natural  resonances 
extracted  (matched  filter  waveform)  from  measured  data  at  bow-on,  stern- 
on  and  abeam  aspects  of  naval  vessels  were  different.  The  stern-on  and 
bow-on  aspect  resonances  were  similar  but  not  the  same  and  both  were 


examination  it  appeared  that  the  aspect  invariance  of  the  resonances 
does  not  hold.  However,  it  was  pointed  out  that  at  these  widely  different 
aspects  different  substructures  of  the  vessels  could  be  excited.  Thus 
for  geometrically  complicated  naval  vessels  al_l_  of  the  resonances  are 
actually  there  but  only  the  observed  subsets  of  these  resonances  are 
excited  at  the  given  aspects.  Similar  effects  can  be  demonstrated  for 
even  such  a simple  target  as  a finite  length  conducting  cylinder.  It 
was  postulated  that  the  observed  complex  natural  resonances  would  actually 
prove  to  be  quasi-invariant,  i.e.,  the  complex  natural  resonances  extracted 
at  the  bow-on  aspect  would  also  permit  vessel  identification  for  some 
range  of  aspects  around  the  bow-on  geometry.  A principal  goal  of  our 
research  during  this  interim  was  to  prove  this  postulate. 

In  the  next  section  of  this  report  measured  multiple  frequency 
scattering  data  are  presented  (representatively)  for  aspects  in  the 
vicinity  of  bow-on,  stern-on  and  abeam  for  two  naval  vessels.  The  next 
section  of  the  report  demonstrates  a clustering  of  the  complex  natural 
resonances  extracted  from  these  data  with  those  resonances  previously 
obtained. 

In  the  next  section  of  the  report  entitled  "Data  Record  Combin- 
ations", the  problem  of  extracting  an  identification-usable  difference 
equation  from  target  records  corresponding  to  various  aspects  is  discussed. 

It  is  shown  that  the  method  used  for  the  data  analysis  in  the  studies, 

Prony's  method  and  clustering,  is  but  one  of  several  different  approaches 
some  of  which  may  offer  more  convenient  and  compact  identification  paramtcers. 
The  identification  results,  i.e,  probabilities  of  target  identification 
are  updated  to  include  the  new  measured  results  and  resultant  "clustered" 
resonances. 


In  the  next  section  results  on  an  encapsulated  source  and  receiver 
system  are  presented.  Using  a worst  case  target  (a  hemisphere  on  a ground 
plane)  data  are  obtained  using  a simple  time  domain  ref lectometry  (TOR) 
system.  The  poles  extracted  from  this  data  are  found  to  agree  quite 
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well  with  those  obtained  by  Stratton  [11]  for  the  perfectly-conducting 
sphere  in  free  space. 

A final  section  of  the  report  summarizes  our  conclusions  to  date 
on  the  feasibility  of  naval  vessel  identification  using  radar  data  and 
outlines  our  intended  research  tasks  during  the  next  interim. 

II.  ASPECT  CLUSTER  DATA  OF  TWO  SHIP  MODELS 

To  test  the  aspect  invariance  postulate  it  was  necessary  to  extract 
resonances  from  model  target  data  for  several  sets  of  aspects  clustered 
about  major  aspects  of  the  models.  The  major  aspects  selected  were  bow, 
broadside  and  stern  incidence  (only  one  broadside  set  of  data  was  obtained 
for  each  ship  since  the  targets  are  symmetric  about  their  centerlines). 

An  angular  increment  of  five  degrees  was  selected  yielding  three  unam- 
biguous views  in  each  direction  from  the  major  aspects.  Thus,  each  target 
was  measured  at  11  aspects  corresponding  to  0°(5°)10o,  80o(5°)100°  and 
170°(5°)180°  (see  Figure  1).  Note  that  the  angle  is  degrees  measured 
from  bow-on,  thus  0°  corresponds  to  bow-on,  90°  to  abeam  and  180°  to 
stern-on  incidence.  It  should  also  be  noted  that  since  a bistatic  radar 
system  was  used,  bow-on  incidence  corresponds  to  the  centerline  of  the 
ship  bisecting  the  antenna  angle.  The  harmonic  frequency  radar  system  and 
ground  plane  range  geometry  are  more  fully  explained  in  references  [1,2]. 

The  range  used  to  obtain  these  data  was  very  similar  to  that  used  previously 
except  for  some  minor  variations  including  placement  of  absorber  panels. 

In  addition,  changes  were  made  in  the  R.F.  system  such  that  more  power 
was  coupled  into  the  higher  order  harmonics*.  For  these  reasons  there 
are  slight  differences  in  the  data  reported  here  and  in  the  previous 
study  Ll]  for  corresponding  aspects. 

*The  major  modification  was  the  use  of  a Varian  ortho-mode  modulator 
in  place  of  the  balanced  modulator  at  X-band  (harmonics  7,  8,  9,  and 
10  of  the  1.085  GHz  fundamental).  The  balanced  modulator  was  used  to 
modulate  the  signals  at  the  20.2  MHz  reference  from  the  network  analyzer 
(Reference  [2],  Appendix  A)  but  did  not  have  the  necessary  bandwidth 
to  cover  all  four  harmonics.  The  use  of  the  orthomode  modulator  has 
significantly  increased  the  output  power  at  harmonics  7,  8 and  9 al- 
though harmonic  10  is  still  relatively  weak. 
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Figure  1.  Radar  interrogation  aspects. 


In  Figures  2a- k and  3a-k  are  given  the  measured  amplitude  data 
for  the  1:500  scale  model  Sverdlov  (a  Russian  Cruiser)  and  the  1:700 
scale  model  Missouri  (a  U.S.  Battleship),  respectively.  Also  shown  on 
those  figures  are  interpolated  spectra  (see  Reference  [1])  which  are 
not  used  at  present.  From  these  data  at  the  harmonic  frequencies,  matched- 
filter  type  responses  [1,2]  were  constructed  as 
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where  6(  6, <J»,p,jnu0)  is  the  square  root  of  the  radar  cross  section  of 


the  target  at  the  n-th  harmonic  as  a function  of  view  angle  (6, 4>)  and 
polarization  (p).  The  principle  advantage  of  this  waveform  (whose  justi- 
fication is  given  in  References  [1]  and  [2])  is  that  only  amplitude  data 
are  required  in  contrast  to  the  ramp- type  responses  used  in  earlier  studies 
which  require  both  amplitude  and  phase.  This  use  of  phase  data  dictates 
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Figure  2j.  1/500  scale  Sverdlov,  175  from  bow-on 
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Figure  3d.  1/700  scale  Missouri,  80°  from  bow-on 
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Figure  4.  1/500  scale  Sverdlov,  0°  from  bow, 

15  poles,  vertical  scale  9.728/1. 


matched  filter  response  that  satisfies  our  present  purpose  it  is  only 
one  of  many  possible  identification  tools.  Further  research  may  point 
tp-otber  more  meaningful  waveforms.  . One  such  .waveform  is  suggested  in 
the  recommendations  section  of  this  report. 
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Using  the  matched-filter  responses,  poles  were  extracted  via  Prony's 
method  [1,2]  using  an  interactive  computer  code  developed  at  this  feci  1 i ty. 
The  poles  were  extracted  from  the  region  indicated  as  1-200  on  Figure 
4 which  corresponds  to  a full  scale  time  span  of  0.230ysec  for  the  1/500 
scale  models  and  0.323psec  for  the  1/700  scale  models.  This  time  span 
is  the  first  half  of  the  period  of  the  fundamental  interrogation  frequency 
and  is  the  largest  unambiguous  portion  of  the  matched  filter  response. 

Since  the  matched  filter  response  results  from  convolving  a step  response 
waveform  with  itself,  this  yields  a symmetrical  waveform  (see  Figure  4, 
solid  line).  In  principle  this  eliminates  the  need  for  phase  data. 

Ideally,  analysis  of  the  complete  waveform  would  yield  frequencies  with 
both  real  and  imaginary  mirror  images,  i.e.,  twice  the  actual  number 
of  frequencies  (poles)  with  half  in  the  right  half  complex  frequency 
plane.  Unfortunately,  the  number  of  poles  needed  in  this  case  (in  the 
neighborhood  of  thirty)  causes  numerical  problems  due  to  the  inversion 
of  a rather  large  matrix.  We  content  ourselves  therefore  with  processing 
only  half  of  the  matched  filter  response  but  recognize  that  use  of  the 
complete  waveform  would  be  desirable.  Note  carefully  that  use  of  the 
full  waveform  should  not  simply  yield  images  of  the  half-waveform  results. 

It  was  not  felt  that  the  required  precision  and  computer  size  could  be 
justified  at  this  point.  We  would  recommend  however  that  future  utilization 
of  our  methods  use  the  complete  matched  filter  waveform.  The  fact  that 
the  waveform  is  symmetric  means  that  the  coefficient  vector  must  also 
be  symmetric.  It  should  be  possible  therefore,  using  a modified  analysis, 
to  utilize  the  complete  matched  filter  response  waveform  without  doubling 
the  size  of  the  coefficient  vector.  It  is  intended  to  study  the  modifi- 
cations necessary  for  this  analysis. 

Since  both  poles  and  residues  are  extracted  with  the  present 
process  LI, 2]  we  are  able  to  reconstruct  the  approximations  to  the  original 
.waveform.  An  example  of  a reconstructed  matched  filter  response  is  given 
by  the  dashed  line  in  Figure  4.  Using  the  pole  extraction  program, 
several  sets  of  poles  were  extracted  from  a given  waveform  by  varying 
the  unspecified  parameters.  Using  each  set  of  poles  and  the  matched 
filter  waveform  a correlation  coefficient  curve  was  calculated  as 


r 


u 


p 

0 


L 


I l-FRM^nAt)  ' FRC^nAt^ 

p"(At)  * 1 - — 


I Fj^nM)  + I Fpr(nA) 


(2) 


RC' 


where  FRM  is  the  matched  filter  response  and  FR£  is  a response  calculated 
using  the  poles,  the  matched  filter  response  and  Corrington's  or  Prony's 
difference  equation  [.1 ,23 - The  pole  set  used  to  specify  the  target  in 
our  library  of  poles  is  defined  to  be  the  one  which  maximizes  the  average 
of  the  correlation  coefficient  over  all  possible  values  of  the  sample 
interval  At  [2].  A summary  of  the  resulting  poles  is  given  in  Tables 
I and  II.  It  should  be  noted  that  the  poles  in  these  tables  have  been 
scaled  such  that  their  imaginary  parts  are  in  mega-radians  per  second. 
Furthermore,  the  poles  are  given  in  terms  of  the  full  scale  frequencies 
for  each  model.  For  a given  harmonic  number  (1  <.  N < 10)  and  a given 
scale  factor  (S  = 500  or  700)  the  radian  frequencies  are 


- 2itN 
“N  ‘ T- 


x 1.085  x 10s 
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where  1.085  GHz  is  the  fundamental  interrogation  frequency.  Since  the 
interrogation  frequencies  should  only  excite  resonances  within  a limited 
region  we  would  expect  the  extracted  poles  to  lie  within  a region  whose 
maximum  is  approximately  given  by  Equation  (3)  with  N equal  to  ten. 

Thus  for 


ufoAX,500  * 136  Mrad/sec 
^X,700  * 97  Mrad/sec* 


(4) 


(5) 


Inspecting  Tables  I and  II  carefully  it  can  be  seen  that  almost  all  the 
poles  having  appreciable  residue  magnitudes  satisfy  these  criteria. 

The  poles  having  imaginary  parts  in  excess  of  the  limits  are  generally 
termed  "curve  fitting  poles"  and  are  typically  associated  with  noisy 
data.  It  is  hoped  that  data  record  combination  techniques  such  as  those 
to  be  discussed  will  help  to  remedy  these  problems  as  an  automated 
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technique  and  also  reduce  the  number  of  entries  in  the  pole  library. 

It  is  stressed  that  ai 1 of  the  parameters  extracted  by  the  program  are 
shown  in  the  Tables.  No  attempt  has  been  made  to  enhance  the  results 
by  removing  obviously  erroneous  poles  nor  do  we  choose  to  plot  only  those 
results  which  obviously  support  our  conjectures.  Admittedly,  averaged 
results  for  the  poles  will  be  enhanced  by  a judicious  choice  of  those 
poles  to  be  included  in  the  averaging  but  Tables  I and  II  are  shown  unaltered 
to  permit  the  reader  to  draw  his  own  conclusions.  We  maintain  that  within 
the  known  limitations  of  the  Prony  method  for  noisy  data  the  extracted 
poles  as  obtained  by  our  approach  do  support  the  quasi-aspect  invariance 
postulate. 

The  next  step  in  the  data  analysis  was  to  treat  the  poles  in  Tables 
I and  II  as  a library  of  entries.  Each  harmonic  data  set  was  then  treated 
as  an  unknown  target  data  set  and  compared  to  all  of  the  pole  library 
entries,  via  the  prediction-correlation  procedure.  The  results  of  this 
test  are  given  in  Table  III.  For  the  harmonic  data  sets  the  prefixes 
"MIS"  and  "SVE"  indicate  the  1/700  scale  Missouri  and  the  1/500  scale 
Sverdlov,  respectively.  The  three  digit  suffix  indicates  the  aspect 
of  the  data  set  in  degrees  from  bow-on  incidence.  For  the  pole  sets 
the  prefixes  "PMI"  and  "PSV"  stand  for  the  best  pole  sets  of  the  Missouri 
and  Sverdlov,  respectively  for  the  indicated  aspect.  The  arrows  indicate 
which  pole  set  gave  the  best  agreement  for  the  criterion  of  maximizing 
the  average  value  of  the  correlation  coefficient,  p",  over  all  possible 
sample  values.  At.  Correct  identifications  are  those  for  which  a data 
set  was  matched  with  poles  extracted  from  the  given  data  set.  A close 
identification  is  one  for  which  the  match  wasn't  with  the  poles  from 
a given  data  set  but  with  a set  of  poles  from  a data  set  within  the 
aspect  group  of  the  data  set  in  question.  An  incorrect  identification 
is  one  in  which  a data  set  was  not  matched  within  the  aspect  group  to 
which  it  belonged.  Thus  even  though  the  aspect  invariance  postulate 
doesn't  appear  to  hold  exactly  it  should  be  possible  to  combine  data 
sets  within  one  of  these  broad  aspect  groups  and  significantly  decrease 
the  number  of  pole  library  entries.  This  prospect  is  enhanced  by  the 
absence  of  mis-identif ications  and  the  similarity  of  some  of  the  pole 
sets  within  the  aspect  groups. 
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TABLE  III 


Harmonic  Data  Set 


Pole  Set 


ABEAM 


STERN 


ABEAM 


STERN 


MIS000 


180 

SVE000 


PMI000  | 

005  >BOW 


090  f ABEAM 


170  j 

175  > STERN 


PSVOOO 


005  >BOW 

010  J 


090  > ABEAM 


175  > STERN 


# Correct 

# Close 

# Wrong 


; 86. 4X 

; 13.6% 
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III.  DATA  RECORD  COMBINATIONS 


For  a given  span  of  real  interrogating  frequencies  where  the  elec- 
trical size  of  the  target  is  not  too  large,  the  target  can  be  characterized 
by  a finite  set  of  complex  natural  resonances  which  are  excitation  invariant. 
Not  all  of  the  resonances  however  are  strongly  excited  at  any  one  aspect 
and  it  may  not  be  possible  to  excite  certain  of  the  resonances  at  that 
aspect.  This  is  particularly  true  for  geometrically  complicated  structures 
such  as  naval  vessels  when  the  resonances  of  Interest  are  physically 
related  to  certain  complex  substructures  of  the  vessel. 

The  details  of  Prony's  method,  which  extracts  a set  of  complex 
natural  resonances  and  residues  from  a given  signal  record,  have  been 
given  [1,2].  With  this  approach  a signal  is  modeled  as  a finite  sum 
of  exponentials* 

N Y„(t-tJ 

f(t)  = l An(e,<|>,p)  e n 0 p(t-t0)  (6) 

and  then  the  complex  natural  resonances 

ir„  = on  i joi„  (7) 

and  the  aspect  and  polarization-dependent  residues 

An(0,4»,p)  = Anr(e,<|>,p)  + j Anq(6,<|»,p)  (8) 

are  extracted  via  a predictor-type  difference  equation  (Prony's  equation) 
[1,2].  There  are  difficulties  with  the  method  in  that  right  half-plane 
poles  (an>0)  and  anomalous  pattern-fitting  poles  can  result.  One  also 
finds  that  the  "poles"  extracted  from  signals  corresponding  to  different 
aspects  of  the  same  target  are  not  precisely  invariant.  Usually  the 
obvious  "pattern  fitting"  poles  are  discarded  and  then  an  averaging  of 
the  remaining  poles  is  used 


♦The  delay  tQ  In  Equation  (6)  is  usually  ignored  mathematically  but 
used  physically. 


40 


(9) 


■J  * 

aa  - Jwa  “ I iJ1  °i  - Jwi 

for  each  obviously  different  pole  location.  The  summand  in  Equation 
(9)  refers  to  various  aspects.  This  approach  works  reasonably  well 
(identification  [3,4])  but  experience  is  helpful  in  deciding  which  pole 
locations  should  be  averaged.  There  are  other  approaches  which  might 
be  attempted.  Suppose  we  have  data  at  I aspects  (i=l,2,...,I)  then  the 
waveforms  could  be  combined  as 


V*> s T ,1,  Vift-V  • <'0) 

Some  preshifting  of  the  waveforms  (T^ ) is  necessary  because  the  usable 
portion  of  each  signal  record  is  generally  dependent  on  aspect.  A pre- 
weighting (d.)  may  also  be  necessary  if  the  amplitudes  of  the  signals 
are  widely  different.  Prony's  method  could  then  be  used  to  extract  one 
set  of  complex  natural  resonances  from  the  averaged  waveform.  Note 
that  averaging  could  reduce  noise  and/or  clutter  problems.  Assuming 
the  form  of  Equation  (6)  for  each  waveform 

1 1 N i „ Y„(t-Ti) 

fa<t)=T.I  di  l A>,*.P)e"  u(t-T,)  . (11) 

i=l  n=l 

Ignoring  the  time-shifting  of  each  waveform  (mathematically) 

f,(t)  =T  I ^ l 02) 

n=l  i=l 

i.e.,  obtaining  an  averaged  waveform  amounts  to  obtaining  a weighted 
averaged  residue  value  for  each  pole  before  processing  by  Prony's  method. 
We  have  not  tested  this  averaged  waveform  approach  as  yet. 

If  Prony's  method  is  applied  individually  to  each  aspect  as  was 
done  here  then  one  obtains 


That  is,  a pole  set  is  obtained  for  each  aspect  and  the  number  of  poles 
in  the  set  can  also  vary  with  aspect  depending  upon  the  criterion  of 
selection.  Note  that  one  also  obtains  the  corresponding  residues 


aJ(0,*,p)  = Ajr(0f*,p)  + j aJj(0,^p)  . (14) 

n=l,2,...N1 

1=1,2,. ..I 

Some,  but  not  all,  of  the  poles  in  Equation  (13)  will  show  the  "clustering" 
referred  to  earlier.  A pole  set  is  then  obtained  from  the  poles  in 
Equation  (13)  using  some  combination  of  machine  and  operator  processing. 

One  reason  for  using  this  approach  for  the  case  of  naval  vessels  is  that 
very  good  success  has  previously  been  obtained  with  "clustered"  poles 
for  underground  targets.  As  is  seen  elsewhere  in  this  report,  poles 
from  the  naval  vessels  do  indeed  show  "clustering",  but  to  date  the 
identification  results  are,  while  not  disappointing,  clearly  not  optimum. 

It  is  not  felt  that  our  processing  is  complete  to  the  extent  that  a 
clustered  pole  approach  should  be  discarded  at  this  point.  We  are,  however, 
beginning  to  look  at  other  methods  for  obtaining  identification  parameters, 
one  of  which  is  the  averaged  waveform  or  averaged  residue  concept  discussed 
earlier. 


The  basic  problem  of  processing  numerous  data  records  from  the 
same  object  is  being  studied  on  a companion  project  [5].  The  main  thrust 
of  these  studies  to  date  is  that  the  most  basic  question  to  be  answered 
is  whether  samples  of  a data  record 

fW<*-n)tt]  - fiJt(N.n)4 

1*1,2, ... I (15)  j 

can  approximately  satisfy  a homogeneous,  linear  difference  equation 

" 1 J 

Jo  ■ s.  , 

m * 0,1,2,...M 

1 = 1,2,. ..I  (16) 
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where  em  is  the  error.  There  are  N+2  different  methods  for  solving  for 
the  coefficients,  an,  [5]  each  of  which  is  a least  squared  error  solution 
(M>2N+1)  and  one  of  which  is  the  well-known  Prony  method.  The  Prony 
method  incidently  does  not  lead  to  the  minimum  totaled  squared  error 
but  this  of  course  is  not  the  only  basis  for  selection  of  a method. 

Equation  (16)  is  slightly  more  general  than  the  equation  given  in  [5] 
in  that  the  increment  on  m is  not  necessarily  the  same  as  on  n.  For 
our  present  purpose,  the  interest  is  in  a suggested  eigenanalysis  method 
for  combining  records  from  several  aspects  by  first  finding  that  difference 
equation  which  best  "fits"  all  of  the  data.  It  has  already  beer,  demon- 
strated that  this  can  be  done  L6,7]  using  chirp-type  short  pulse  response 
waveforms  from  aircraft  targets.  It  is  intended  to  test  the  naval  vessel 
data  using  this  approach. 

From  a pole  or  complex  natural  resonance  viewpoint  the  identification 
problem  is  complicated  by  the  fact  that  we  know  from  present  results 
that  for  a given  vessel  at  least  three  distinct  pole  subsets  can  be  defined 
corresponding  to  aspects  in  the  vicinity  of  bow-on,  stern-on  and  abeam. 

It  is  not  yet  clear  if  these  pole  subsets  should  be  treated  as  three 
separate  subvessels,  all  of  which  identify  che  same  vessel,  or  if  the 
subsets  should  be  combined  to  characterize  the  vessel  and  the  processing 
permitted  to  reject  the  unused  or  unexcited  poles.  To  date  only  the 
subvessel  approach  has  been  tested  and  then  only  incompletely.  Given 
that  a set  of  complex  natural  resonances  can  be  used  to  characterize 
a given  vessel  (we  absorb  the  as  yet  unanswered  question  of  several  sub- 
sets vs.  one  larger  set)  it  follows  that  there  are  many  possible  ways 
to  identify  the  target.  The  prediction-correlation  identification  pro- 
cess developed  here  has  worked  better  than  any  other  method  we  have  tried. 
It  has  the  additional  very  real  advantage  that  it  is  not  necessary  to 
extract  poles  from  some  unknown  target  record.  This  does  not  mean  of 
course  that  a better  method  of  identification  processing  might  not  be 
found  in  the  future. 
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IV.  THE  ENCAPSULATED  SOURCE  AND  RECEIVER 


A major  problem  with  the  identification  techniques  presented  is 
obtaining  characteristic  waveforms  for  the  targets  in  question.  Once 
this  is  accomplished,  then  using  Prony's  method  or  some  similar  method 
the  poles  may  be  extracted  and  identification  can  be  accomplished.  For 
this  project  the  signaling  waveforms  have  been  produced  via  Fourier 
synthesis  of  discrete  frequency  measurements  of  scaled  model  targets. 

As  we  attempt  to  go  to  higher  frequencies  practical  limits  for  an  experi- 
mental system  are  quickly  reached.  As  explained  earlier  we  have  modeled 
the  ocean  environment  by  a highly  conducting  ground  plane.  Even  at  the 
frequencies  being  used  now  this  model  is  in  question  and  the  problem 
will  become  more  serious  at  higher  frequencies.  Also,  as  larger  models 
are  used  the  problem  of  setting  up  a measurement  system  which  can  accurately 
measure  the  scattered  field  from  the  targets  becomes  increasingly  difficult. 

With  larger  targets  the  problems  of  near  field  effects  become  more 
serious.  The  required  distances  to  alleviate  this  problem  would  make 
an  experimental  system  impractical.  This  problem  would  become  even  more 
formidable  if  it  becomes  necessary  to  look  at  the  targets  from  many  different 
aspects.  It  has  been  suggested  [9,10]  that  a measurement  system  utilizing 
an  encapsulated  source  and  receiver  could  be  used  to  obtain  characteristic 
waveforms  of  the  targets.  Obviously  encapsulating  the  source  and  receiver 
does  not  correct  the  sea  model  (ground  plane)  at  high  frequencies.  The 
point  however  is  that  encapsulated  source-receiver  techniques,  when  fully 
developed,  could  be  applied  to  full  scale  targets  in  situ.  In  Figure 
5 is  shown  a simplified  diagram  of  a system  using  this  concept.  This 
system  is  similar  to  those  developed  for  use  in  the  area  of  underground 
radar  [12].  A short  pulse  source  is  used  to  excite  the  target  and  a 
sampling  scope  is  used  to  record  the  resulting  forced  response  and  transient. 
Typically  the  forced  response  is  many  times  greater  in  magnitude  than 
the  transient  portion.  Thus  if  sensitivity  in  the  transient  region  is 
not  to  be  sacrificed  then  some  method  must  be  found  to  protect  the  sampling 
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Figure  5. 


head.  In  the  underground  radar  system  this  is  done  by  using  crossed 
dipoles  as  the  excitation  and  receive  antennas.  The  receiver  dipole, 
being  significantly  decoupled  from  the  excitation  dipole,  picks  up  a 
transient  response  which  is  large  enough  to  be  accurately  read  but  rejects 
the  major  portion  of  the  target's  forced  response.  In  the  case  of  the 
encapsulated  source  and  receiver  it  would  be  possible  to  use  two  antennas 
although  one  antenna  would  be  more  convenient  for  a full  scale  system. 

However  the  crossed  dipole  system  or  an  analogous  type  of  isolating  system 
would  be  difficult  to  implement.  This  then  complicates  the  measurement 
system  which  is  contrary  to  our  purpose. 

In  Figure  6 is  shown  a simplified  system  we  have  tried  which  utilizes 
a modified  version  of  the  encapsulated  source  and  receiver  concept. 

This  system  uses  a time-domain-ref lectometry  (TDR)  unit  which  is  connected 
via  a six  foot  section  of  coaxial  cable  to  a stub  antenna  on  a target. 

The  TDR  unit  sends  out  a series  of  short  pulses  (in  this  case  a tektronix 
TDR  unit  has  been  used  with  a pulse  length  of  approximately  120  picoseconds). 
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Figure  6. 

A sampling  scope  is  used  to  sample  the  reflected  energy  from  down  the 
line.  The  output  is  a continuous  waveform  showing  reflection  coefficient 
as  a function  of  time  (or  range  if  the  velocity  of  propagation  in  the 
medium  is  known).  In  the  case  of  the  hemisphere  on  the  finite  ground 
plane  shown  in  Figure  6 we  are  interested  in  the  region  of  the  response 
due  to  the  stub-hemisphere  connection  and  beyond.  In  this  broad  region 
we  can  expect  several  subregions  of  interest.  First  there  will  be  a 
sharply  increasing  region  due  to  the  transition  from  the  cable  center 
conductor  surrounded  by  a dielectric  to  the  stub  antenna  which  is  in 
free  space  (see  Figure  7).  Next  there  will  be  a region  of  special  interest. 
This  should  be  a series  of  oscillations  in  impedance  due  to  reflections 
from  the  ground  plane.  Note  that  the  period  of  the  oscillation  is  dependent 
on  the  size  of  the  sphere  and  is  therefore  characteristic  of  It.  We 
will  also  expect  a return  from  the  edge  of  the  ground  plane  which  should 
be  far  enough  away  to  be  Isolated  from  the  hemisphere  oscillations. 
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Figure  7.  Example  of  TDR  response. 

There  will  also  be  oscillations  due  to  returns  from  the  end  of  the  stub 
antenna.  If  the  antenna  is  too  short  then  these  oscillations  could  inter- 
fere with  those  due  to  the  hemisphere.  Lastly  there  will  be  a constant 
section  which  is  indicative  of  complete  transition  to  free  space. 

In  Figure  8a,b,c  are  shown  the  TDR  responses  of  a 7.05  inch  diameter 
hemisphere  on  a 48  x 48  inch  square  metal  plate.  The  stub  lengths  are 
0.78  inches,  3.00  inches  and  29.25  inches  respectively.  In  Figures  8a 
and  8b  it  is  easily  seen  that  the  oscillations  due  to  the  antennas  tend 
to  cover  up  those  of  the  hemisphere.  In  Figure  8c  the  response  of  the 
hemisphere  is  easily  distinguishable  from  the  returns  due  to  the  edge 
of  the  ground  plane.  For  this  reason  all  processing  has  been  done  on 
the  response  in  Figure  8c.  This  method  of  isolating  the  target  response 
from  the  antenna  response  is  easily  done  for  the  case  of  the  small  hemisphere 
but  whould  be  impractical  for  large  targets.  It  is  felt  that  it  might 
be  possible  to  subtract  this  out  using  a difference  equation  approach 
and  careful  modeling  of  the  cable-antenna  transition.  Further  work  is 
anticipated  in  this  area. 
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DATA  ID  : HEMISPHERE,  2.9'STUB,  H-l^-78 
6*  CABLE. 10  NS, AVE=30 
NSTART  : 73 

NLAST  : 202 

* POLES  : 15 


Figure  8c. 

Included  with  Figure  8c  is  a waveform  constructed  from  poles  ex- 
tracted via  Prony's  method.  In  this  case  the  fit  to  the  waveform  was 
very  good  and  over  the  interval  of  application  the  measured  and  recon- 
structed waveforms  are  virtually  identical  (the  reconstructed  waveform 
is  given  by  the  dashed  line).  By  varying  the  number  of  poles  requested 
and  the  interval  ([2],  Section  V)  many  pole  pairs  were  extracted.  The 
results  of  this  operation  are  shown  as  a plot.  Figure  9,  of  the  pole 
positions  in  the  complex  plane.  It  should  be  noted  that  there  is  a 
definite  clustering  effect  of  the  poles.  Also  there  are  many  spurious 
points  that  can  be  associated  with  curve  fitting  poles  which  generally 
have  very  small  residues  associated  with  them.  Also  plotted  in  Figure 
9 are  the  pole  positions  for  the  7.05  inch  sphere  as  obtained  by  Stratton 


Figure  9.  Poles  of  a 7.05"  hemisphere  via  Prony's  method. 

[11].  Due  to  the  stub  excitation  we  would  expect  to  excite  the  transverse 
magnetic  (TM)  resonances  of  the  spherical  body.  The  effect  of  the  ground 
plane  will  be  to  short  out  the  even  mode  resonances.  Thus  we  will  expect 
to  see  the  odd,  TM  modes.  From  Figure  9 it  can  be  seen  that  there  is 
a definite  clustering  of  the  poles  about  the  TM1  and  TM^  resonances. 

There  is  quite  a spread  in  the  real  parts  of  the  poles  but  fairly  good 
agreement  for  the  imaginary  parts.  For  the  TMg  pole  there  is  no  cluster 
of  poles  near,  however,  there  is  a cluster  at  a slightly  greater  frequency 
with  a real  part  approximately  40%  of  that  of  the  true  pole.  This  spread 
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of  the  real  parts  has  been  seen  before  when  Prony's  method  has  been  used 
on  noisy  data.  It  is  felt  that  these  results  are  very  promising  due 
to  the  nature  of  the  sphere  which  is  almost  a worst-case  target.  This 
is  due  to  the  high  damping  of  the  poles  which  results  in  very  few  oscil- 
lations compared  to  targets  such  as  thin  wires  which  show  appreciable 
ringing  for  several  cycles.  For  this  reason  it  is  planned  to  obtain 
measurements  on  other  targets  such  as  the  prolate  spheroid.  This  will 
also  be  cut  in  half  and  measured  on  the  ground  plane.  Due  to  its  long 
slender  shape  hopefully  this  target  will  ring  more  and  be  better  suited 
to  pole  extraction  techniques  like  Prony's  method. 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  can  be  drawn  from  the  factual  data  sections 
of  this  report. 

1.  The  feasibility  of  naval  vessel  identification  via  complex 
natural  resonances  associated  with  substructure  features  of 
the  vessel  has  been  demonstrated. 


Identification  capability  has  been  extended  to  arbitrary  aspects 
in  the  vicinity  of  bow-on,  stern-on  and  abeam,  i.e.,  the  quasi- 
invariant postulates  regarding  previous  results  have  been 
proven. 

By  implication,  substructure  complex  natural  resonances  and 
prediction-correlation  processing  can  be  extended  to  arbitrary 
aspect  angles  but  the  number  of  resonance  aspect  positions 
for  a given  vessel,  while  small,  is  unknown. 
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4.  Interactive  programs  based  on  Prony's  method  have  been  developed 
for  extracting  poles  and  residues  from  noisy  transient  response 
signals.  The  programs  automate  the  selection  of  sample  interval 
and  pole  number.  These  programs,  while  adequate  for  our  present 
task,  still  display  many  of  the  known  problems  associated 
with  Prony's  method. 

5.  Other  methods  for  extracting  a set  of  poles  and  residues  from 
several  transient  records  corresponding  to  different  target 
aspects  have  been  suggested.  Certain  of  these  are  being  pursued 
on  companion  (JSEP)  contracts. 

6.  Ten  frequency  harmonic  scattering  data  on  two  naval  vessel 
classes  (cruiser  and  battleship)  have  been  obtained  for  a 
variety  of  aspects  spanning  the  bow-on,  stern-on  and  abeam 
positions.  These  data  served  as  test  inputs  for  the  extended 
identification  procedures. 

7.  For  two  naval  vessel  classes  and  aspects  in  the  vicinity 
of  bow-on,  stern-on  or  abeam  probabilities  of  correct  identi- 
fication of  86  percent  and  probabilities  of  incorrect  identifi- 
cation of  14  percent  have  now  been  obtained.* 

1 

8.  The  feasibility  of  using  encapsulated  source  and  receiver 
methods  to  obtain  natural  resonance  rich  transient  signals 
has  been  demonstrated.  These  methods,  when  refined,  could 
be  used  to  obtain  the  complex  natural  resonances  of  large 
geometrically  complicated  structures  in  situ.  Test  results 
for  conducting  hemispherical  geometries  on  a ground  plane 
were  extremely  promising,  i.e.,  the  measured  resonances  were 
in  good  agreement  with  the  known  resonances. 


♦These  are  conservative  estimates  of  identification  probabilities. 
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A fair  summary  of  our  present  progress  with  regard  to  the  radar 
identification  of  naval  vessels  can  be  stated  in  terms  of  what  appears 
to  be  possible  now  with  sophisticated  full  scale  radar  systems.  Given 
a radar  system  capable  of  harmonically  sampling  the  radar  cross  section 
of  a naval  vessel  in  situ  over  a frequency  decade  with  the  highest  fre- 
quency less  than  25  MHz,  correct  identification  of  the  vessel  with  prob- 
abilities in  the  range  70  to  80  percent  are  postulated.  While  these 
correct  identification  probabilities  are  low,  they  are  tempered  by  the 
fact  that  the  probabilities  of  incorrect  identification  are  remarkably 
low.  In  principle,  these  results  are  independent  of  the  aspect  and  elevation 
angles  of  the  radar  and  the  target  range.  It  is  assumed  that  a vertically 
polarized  radar  is  used.  The  frequency  range  of  the  system  is  such  that 
the  sea  state,  minor  structural  alterations  of  the  vessel  and  the  motion 
of  the  vessel  are  not  factors  which  will  contaminate  or  distort  the 
identification  process.  A major  consideration  js  obtaining  the  needed 
sampling  of  transient  response  waveforms  of  the  vessel  from  which  the 
identification  parameters  (substructure  complex  natural  resonances)  can 
be  extracted.  At  the  moment  identification  parameters  on  only  two  vessels 
(see  Section  V)  are  available  and  these  are  incomplete.  Note  carefully 
that  transient  data  for  extraction  of  the  identification  parameters  need 
not  be  at  any  particular  range,  aspect  or  elevation,  but  must  span  dis- 
cretely all  possible  major  aspects.  At  the  very  least  this  would  include 
the  bow-on,  abeam  and  stern-on  regions  plus  angles  45°  and  135°  from 
bow-on.  These  conjectures  are  based  on  bistatic  (10°)  measurements  but 
no  significant  alterations  for  true  monostatic  systems  are  anticipated. 

Perhaps  the  most  direct  and  immediate  need  for  this,  as  well  as 
other,  radar  identification  programs  is  multiple  frequency  scattering 
data  over  appropriate  spectral  regions  (very  roughly  hull  lengths  of 
1.5  to  10  or  15  wavelengths).  Some  care  has  been  taken  to  ensure  that 
our  scattering  models  are  sound  both  geometrically  and  electromagnetically 
[1]  and  that  the  required  frequencies  for  identification  are  at  least 
feasible  if  not  optimum  from  a practicality  viewpoint.  Unfortunately 
it  is  precisely  this  spectral  range  which  dictates  experimental  data 
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as  neither  computer  models  nor  asymptotic  theories  can  be  used  in  their 
present  forms.  Multiple  frequency  measurements  of  this  type  are  tedious 
and  time  consuming.  Additional  details  (frequencies)  are  also  needed 
if  the  substructure  features  associated  with  a particular  identification 
parameter  are  to  be  isolated.  In  essence,  simple  models  and  measurements 
may  be  adequate  for  feasibility  studies  but  we  are  beyond  that  stage. 

It  is  desired  to  test  our  identification  procedures  over  as  broad  a data 
base  as  possible  (spanning  at  least  several  vessel  classes  for  a variety 
of  aspects,  elevations  and  ranges)  while  at  the  same  time  refining  our 
processing  procedures.  The  recommendations  listed  below  are  directed 
toward  these  goals. 

1.  A continual  testing  of  the  available  data  base  (two  vessels) 
using  different  criteria  to  select  the  size  and  content  of 
the  identification  parameter  library.  These  procedures  can 

and  will  be  expanded  to  include  new  data  as  they  become  available. 

2.  New  and  unique  discrete-swept  scattering  measurements  of  the 
naval  vessels  in  situ  now  available.  These  are  essentially 
swept  frequency  measurements  but  the  frequency  is  stabilized 
and  the  background  corrected  at  many  discrete  frequencies 
across  the  frequency  span.  For  example,  in  the  range  2.2 

to  4.0  GHz  (model  frequencies)  some  201  spot  measurements 
are  made.  These  measurement  procedures  have  been  automated 
and  tested  for  free  space  targets  and  are  presently  being 
adjusted  to  the  naval  vessel  range. 

3.  From  the  measured  data  in  2)  above,  synthetic  pulse  response 
waveforms  can  be  synthesized.  We  propose  to  obtain  resolutions 
such  that  the  significant  scatterers  of  the  vessel  can  be 
isolated  and  the  dominant  substructure  controlling  resonances 
identified.  The  pulse  response  waveform  envelopes  will  also 
permit  certain  nonphysical  but  possibly  useful  identification 
tests  to  be  made.  Details  of  the  discrete-swept  reflectivity 
facility  are  given  in  References  [6,7],  as  are  examples  of 
envelope  processing. 
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Analytical  models  for  calculating  the  radar  return  from  in 
situ  vessel  models  are  beginning  to  emerge  [8].  While  the 
ship  models  are  still  relatively  simple  and  the  ship  scattering 
theory  confined  to  first  order  diffraction,  a nonsmooth  sea 
surface  and  ship-sea  iteraction  are  included.  Thus  the  basic 
model  would  appear  to  be  a correct  basis  for  more  sophisticated 
models.  We  would  propose  a beginning  look  at  these  models, 
perhaps  by  comparing  our  existing  measurements  with  some  simple 
computations.  A carefully  developed  model  which  realistically 
(interactions)  includes  a model  of  a nonsmooth  sea  is  a 
desirable  and  economically-wise  alternative  to  very  expensive 
and  time-consuming  measurements  which  cannot  evaluate  changing 
sea  states.  In  addition,  such  models  when  developed  could 
provide  a significant  data  base.  This  item  is  viewed  as  a 
long  term  multi-year  objective. 

Certain  moderate  scale  encapsulated  source  and  receiver  experi- 
ments should  be  made.  Large  models  (dimensions  of  a meter 
roughly)  of  simple  models  such  as  cubes  or  paral lelpipeds 
can  be  constructed  with  reasonable  cost.  The  use  of  such 
models  would  permit  bistatic  arrangements  of  source  and  receiver 
probes  thereby  removing  the  TOR  restrictions  on  the  shock- 
type  excitations.  Also,  the  use  of  screened  rooms  for  the 
measurements  should  not  be  necessary.  One  or  two  experiments 
of  this  type  are  recommended  knowing  in  advance  that  some 
probe  development  work  will  be  necessary  and  that  the  complex 
natural  resonances  of  such  targets  are  not  known  analytically. 

Current  studies  on  companion  contract  (project  710816)  show 
promise  of  developing  resonance  extracting  procedures  which 
are  less  sensitive  to  noisy  transient  signals.  Additional 
study  along  these  lines  is  only  recommended  as  found  necessary 
to  adjust  general  techniques  to  the  particular  data  on  this 
contract. 
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